February 1, 2008 4:38 Proceedings Trim Size: 9in x 6in 



sewm04'ejiri 



o 
o 

(N . NUMERICAL STUDY OF THE EQUATION OF STATE FOR 

bJO 1 TWO FLAVOR QCD AT NON-ZERO BARYON DENSITY 

< 



S. EJIRJf C.R. ALLTON^ M. DORINGf S.J. HANDS? 3 O. KACZMAREKf 
F. KARSCHf E. LAERMANNf and K. REDLICH ac 

a Fakultat fur Physik, Universitdt Bielefeld, D-33615 Bielefeld, Germany 
^ ' b Department of Physics, University of Wales Swansea, Singleton Park, 

\Q ', Swansea, SA2 8PP, U.K. 

^Institute of Theoretical Physics, University of Wroclaw, PL-50204 Wroclaw, 
' Poland 

00 " 

o 

i We discuss the equation of state (EoS) for two flavor QCD at non-zero temperature 

and density. Derivatives of In Z with respect to quark chemical potential fi q are 
^ ' calculated up to sixth order. From this Taylor series, the pressure, quark number 

density and associated susceptibilities are estimated as functions of temperature 
and fig. It is found that the fluctuations in the quark number density increase 
in the vicinity of the phase transition temperature and the susceptibilities start 
to develop a pronounced peak as fj, q is increased. This suggests the presence of 
a critical endpoint in the (fi q ,T) plane. Moreover, we comment on the hadron 
resonance gas model, which explains well our simulation results below T c . 



X 



H 

CS ' 1. Introduction 



It is important to study QCD at high temperature and density by numerical 
simulations of lattice QCD. In particular, studies of the equation of state 
(EoS) can provide basic input for the analysis of the experimental signatures 
for quark-gluon plasma formation. We study the EoS in the low but non- 
zero baryon number density regime, interesting for heavy-ion collisions. 
The simulation of QCD at non-zero fi q is known to be difficult. However, 
studies based on a Taylor expansion with respect to fi q turned out to be an 
efficient technique to investigate the low density regime ' . 

Thermal fluctuations near the critical temperature T c are important in- 
formation for the understanding of the QCD phase diagram. Baryon num- 
ber fluctuations are expected to diverge at the endpoint of the first order 



'Presented by S. Ejiri. This work is supported by BMBF grant No.06BI102, DFG grant 
KA 1198/6-4, PPARC grant PPA/a/s/1999/00026 and KBN grant 2P03 (06925). 
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chiral phase transition line. The electric charge fluctuation in heavy-ion col- 
lisions is one of the most promising experimental observables to identify the 
critical endpoint. The fluctuations of quark number, isospin number and 
electric charge are estimated by the corresponding susceptibilities, Xq,Xi 
and XCi which are given by \q — d 2 p/d/j, 2 and \i ~ d 2 p/dfj,j 1 where 
fi q = (p u + /id)/2, /i/ = (nu — A i d)/2, and fJ, u [d] is the chemical potential for 
the u[d\ quark. For fi u = [14, the charge susceptibility is \c — Xg/36+x//4. 

In this study, we investigate the pressure p, quark number suscep- 
tibility Xq an d isovector susceptibility \i by calculating the Taylor ex- 
pansion coefficients of p, c n and c^, up to sixth order in /i, which are 
defined by p/T 4 = (lnZ)/(VT 3 ) = En=o c n(Vq/T) n and Xq [xi]/T 2 = 
En=2 n ( n - ^Cn[ci}(^ q /T) n - 2 for M „ - fx d . where c„ = (d n \nZ/d^)/n\ 
and = (d n \nZ/dfiq~ 2 dfi 2 )/n\. Moreover, we compare these results with 
the prediction from the hadron resonance gas model 3 ' 4 . 

2. Quark gluon gas or hadron resonance gas 

We expect the equation of state approaches that of a free quark-gluon gas 
(Stefan-Boltzmann (SB) gas) in the high temperature limit. The coefficients 
in the SB limit for Nf = 2, fij = are well known as C2 = = 1, C4 = c| = 
1/(2tt 2 ) and c„ = c T n = for n > 6. 

On the other hand, in the low temperature phase QCD is well modelled 
by a hadron resonance gas. If the interaction between these hadrons can 
be neglected, the pressure is obtained by summing over the contributions 
from all resonance states of hadrons. The contribution to p/T 4 from a 
particle which has mass m^, baryon number Bi and third component of 
isospin number I 3i is given by 

where K2 is the modified Bessel function, and 77 is 1 for mesons and — 1 
for baryons. Moreover, since rrii/T 3> 1 for all baryons and K2(x) « 
\J n/2x exp(— x), the first term of I is dominant in the baryon sector. There- 
fore, the pressure can be written by p/T 4 = G(T) + F(T) cosh(3/Xg/T) 
for /ij = 0, where G(T) and F(T) are the mesonic and baryonic compo- 
nents of p/T 4 at (i q = 0, respectively 4 . Similarly, we obtain Xq/T 2 = 
9F(T)cosh(^q/T) and X i/T 2 = G 1 (T) + F 1 (T) cosh(3(i q /T). Here, the 
mesonic component for Xq is zero because mesons (Bi = 0) are independent 
of /i q . Therefore, C4/C2 = 3/4, cq/c^ = c\/c\ = 3/10 in the region where 
the non-interacting hadron resonance gas provides a good approximation. 
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Figure 1. The ratio of Taylor expansion coefficients for \ q (left) and \I (right). 



We investigate these coefficients. Simulations arc performed on a 16 3 x 4 
lattice for 2 flavor QCD. The pA improved action 5 is employed at ma = 0.1, 
which gives mps/mv = 0.7. The number of configurations is 1000-5000 
for each T. We use the random noise method 1,2 . The results for c n+ 2/c n 
and cf i+2 /c^ are shown in Fig. 1. We find these results are consistent with 
the prediction from the hadron resonance gas model for T/T c < 0.96 and 
approach the SB values, i.e. C4/C2 = cl/c^ = l/2ir 2 , cq/c4 = Cg/c| = 
0, in the high temperature limit. These results suggest that the models 
of free quark-gluon gas and hadron gas seem to explain the behavior of 
thermodynamical quantities well except in the narrow regime near T c . 

Moreover, if the data is consistent with the hadron resonance gas model, 
we can separate \i m to the contributions from isotriplct mesons and 
baryons, since G 1 + F 1 = 2c\ and 9F 1 /2 = I2C4. Isosinglet mesons do 
not contribute to Xi because 1^ — 0. Also, 9F = 2c2 from the equation of 
Xq- Furthermore, for the case of m, 3> T, the first term (I = 1) in Eq.(l) 
is dominant. This approximation holds for our present simulations where 
m^/T ~ 4-5 but will no longer be valid with physical quark mass values. 
For our present simulation parameters, the ratio of 9 2 (p m i/T 4 )/9(/z//T) 2 
and Pmi/T 4 equals (2/ 3i ) 2 for I 3i = { — 1,0, 1}. We thus find for the con- 
tribution of the triplet meson part: G — 3G 1 /8. We then can estimate the 
contribution of the isosinglet mesons to p/T A by comparing F + G with the 
total pressure. In Fig. 2 (left), we plot the baryon component F = Ic-ijQ 
(square), the triplet meson component G — Sc^/A ~ c\ (circle) and the 
total (diamond). The solid line is the pressure obtained with the integral 
method in Ref. 6. The dashed lines in Fig. 2 (left) represent the predictions 
from the hadron resonance gas model for the mass parameter of our simula- 
tions. The resonance states for this quark mass are adjusted by the method 
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Figure 2. Pressure at fi q = as a function of T (left), and difference of pressure at 
non-zero p. q /T from /i q = for each fixed (J-q/T (right). To is T c at /j, q = 0. 

described in Ref. 3. The simulation results for p/T 4 obtained with the inte- 
gral method, together with the triplet meson and baryon contributions, are 
surprisingly well reproduced by this model calculation. The result of the 
singlet meson contribution may explain the difference in the pressure ob- 
tained from the integral method and that from the Taylor expansion which 
neglects the singlet part. 

3. Pressure and susceptibilities at fi q ^ 

Next, we calculate pressure and quark number susceptibility in a range of 
< HqjT < 1, using the data of c„ and c£; A(p/T 4 ) = p{T 7 fi q )/T 4 - 
p(T,0)/T 4 = cMTf + c^ q /T) 4 + c 6 ( N /Tf, X ,/T 2 = 2c 2 + 
I2ci(p, q /T) 2 + 30cs(iJ, q /T) 4 , and the corresponding equation for \i- 

We draw A(p/T 4 ) for each fixed [i q jT in Fig. 2 (right) and find that 
the difference from p| M?= o is very small in the interesting regime for heavy- 
ion collisions, \i q jT rj 0.1 (RHIC) and fJ. q /T « 0.5 (SPS), in comparison 
with the value at /i q = 0, e.g. the SB value for 2 flavor QCD at /i g = 0: 
p SB /T 4 ~ 4.06. The effect of a non-zero quark density on the pressure at 
H q /T = 0.1 is only 1%. Comparing with previous result up to 0(fj, 4 ) 2 , 
the effect from 0(/j, q ) seems to be small in the range we investigated. Also, 
the result is qualitatively consistent with that of Ref. 7 obtained by the 
reweighting method. 

The data connected by solid lines in Fig. 3 are the results for Xq (left) 
and xi (right) obtained by this method up to 0(/i 4 ). We also calculated 
the susceptibilities by another two methods. Dot-dashed lines are from the 
hadron resonance gas model, using F, F 1 and G 1 we obtained. This turned 
out to be a good approximation for T/T c < 0.96. Dashed lines are calculated 
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Figure 3. Quark number (left) and isovector (right) susceptibilities as a function of T 
for each fixed /J-q/T. To is T c at fj, q = 0. 



by the reweighting method using an approximation discussed in Ref. 1. Here 
the truncation error is 0(n®) but the effect from higher order terms of \i q is 
partially included. The sign problem in the reweighting method is serious 
when the complex phase fluctuations of the quark determinant are large at 
large fJ- q /T. Hence we omitted data at which the standard deviation of the 
complex phase is larger than tt/2. The difference among the three results 
is caused by the approximation in the higher order terms of n q /T. 

Since the statistical error of cq is still large near T c , the peak of Xq 
is not clear. However, as seen in Fig. 1, cq changes its sign at T c . This 
means the peak position of Xq moves left, which is corresponding to the 
change of T c as a function of [i q . T c ([i q /T — 1)/T c (fi q /T = 0) in Ref. 1 
is about 0.93. Moreover, the strong enhancement of \q with increasing fi q 
near T c , obtained in Ref. 2, is observed also in the analysis up to 0(fi q ). 
This suggests the presence of a critical endpoint in the (T, (j, q ) plane. On 
the other hand, xi i n Fig. 3 does not show any singular behavior. This is 
consistent with the sigma model prediction that only isosinglet degrees of 
freedom become massless at the critical endpoint 8 . 
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